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Abstract 

We apply our recent results on orbital dynamics around a mass- 
varying central body to the phenomenon of accretion of Dark Matter- 
^ . assumed not self-annihilating-on the Sun and the major bodies of the 

qh| solar system due to its motion throughout the Milky Way halo. We 

^ . inspect its consequences on the orbits of the planets and their satellites 

over timescales of the order of the age of the solar system. It turns 
out that a solar Dark Matter accretion rate of w 10^^^ yr^^, inferred 
from the upper hmit AM/M = 0.02 - 0.05 on the Sun's Dark Matter 
. content, assumed somehow accumulated during last 4.5 Gyr, would 

' have displaced the planets faraway by about 10"'^ — 10^ au 4.5 Gyr 

On ■ ago. Another consequence is that the semimajor axis of the Earth's 

orbit, approximately equal to the Astronomical Unit, would undergo 
a secular increase of 0.02 — 0.05 m yr~^, in agreement with the latest 
observational determinations of the Astronomical Unit secular increase 
. of 0.07 ±0.02 m yr~^ and 0.05 m yr^^. By assuming that the Sun will 

' continue to accrete Dark Matter in the next billions year at the same 

rate as in the past, the orbits of its planets will shrink by about 10~^ — 
10^ au (w 0.2 — 0.5 au for the Earth), with consequences for their fate, 
. especially of the inner planets. On the other hand, lunar and planetary 

' ephemerides set upper bounds on the secular variation of the Sun's 

gravitational parameter GM which are one order of magnitude smaller 
than « 10^^^ yr~^. Dark Matter accretion on planets has, instead, 
less relevant consequences for their satellites. Indeed, 4.5 Gyr ago 
their orbits would have been just 10^^ — 10^ km wider than now. Dark 
Matter accretion is not able to explain the observed accelerations of the 
orbits of some of the Galilean satellites of Jupiter, the secular decrease 
of the semimajor axis of the Earth's artificial satellite LAGEOS and 
the secular increase of the Moon's orbit eccentricity. 
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1 Introduction 

Dark Matter (DM) was postulated long ago to explain the discrepancy be- 
tween the observed kinematics of some components of astrophysical systems 
like clusters of galaxies [Ij and spiral galaxies [21 13] , and the predicted one 
on the basis of the Newtonian dynamics and the matter directly detected 
from the emitted electromagnetic radiation (visible stars and gas clouds). 
Zwicky [1] postulated the existence of undetected, baryonic matter; today 
[5] it is believed that the hidden mass is constituted by non-baryonic, weakly 
interacting particles in order to cope with certain issues pertaining galaxy 
and galaxy clusters formation [6] and primordial nucleosynthesis [7]. On 
cosmological scales, DM accounts for about 23% of the mass-energy content 
of the Universe [8]. 

It is important to perform tests of the DM existence in other scenarios, 
independently of the phenomena themselves for which its existence was orig- 
inally postulated. Otherwise, it could not get rid of the unpleasant status 
of ad-hoc hypothesis. 

Our solar system can be regarded as one of such different laboratories. 
The first arguments in support of invisible (standard) matter in solar system 
were presented in 1932 by Oort [9j in connection with the dynamical density 
in the vicinity of the Sun's Galactic orbit. Oort concluded that the total 
mass of nebulous or meteoric matter near the Sun is less than 3 x 10"'^^ g 
cm~^. Nowadays, it is phenomenologically evaluateci3 from that the upper 
bound on local DM density in the solar system 



is several orders of magnitude higher than the mean Galactic value [10] 



so that it could be imagined that the solar system is surrounded by a local 
subhalo, a suggestion reinforced by recent simulations |Tl]. Several processes 
have been postulated to clump DM in the solar system [HI [I3l [Ml [I5l [I6l 
[T71 [T8l [19] : according to Ref. [Sni 121], the existence of the solar system 
itself might be evidence for a local subhalo. However, it must be noted that, 
according to Ref. [22] , densities of DM as large as that of eq. (Jlj) could not be 
reached because of the inverse process of capturing DM from the Milky Way 




(1) 




^See later. 
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halo. According to it, whenever the density of DM bound to the solar system 
gets big, the process of throwing WIMPs out of it will be in equilibrium 
with capture from the Galactic halo. This fact has been proven by using 
the Liouville's theorem, which says that the phase space density (after a 
long time, when equilibrium has been reached) of the solar system bound 
WIMPs is the same as in the halo \22\ [T5] . Such arguments have recently 
been reviewed in [23]. Several studies have been dedicated to placing bounds 
on the local distribution of DM in our solar system from orbital motions of 
natural major and minor bodies and artificial probes [201 EH ESI ESI EZl EH 
ESlSOlEIlESllMllSlEaEQlIMlEZj; actually, the upper bound of eq. (H} 
comes just from such investigations. Propagation of electromagnetic waves 
has been used as well |38|. It has also been suggested that some proposed 
and approved space-based missions for fundamental physics may be used to 
search for local distributions of DM as well \39\ HOl E] . Apart from directly 
affecting orbital motions and electromagnetic waves propagation, DM may 
also have other effects on the solar system's bodies: for example, it may 
concur to form main-sequence stars like the Sun [32] and accrete on the 
major planets [431 121j. and affect their internal heats [441 145j . According 
to Ref. [43] ■ during its journey along the Galaxy since its birth 4.5 Gyr ago 
the solar system would have encountered about 2O3M0 of DM. In Ref. ^42| 
it is stated that the upper limit to the amount of gravitating DM in the 
Sun would be 2 — 5% of the total solar mass; by assuming that it somehow 
accumulated during the Sun's lifetime, it is possible to infer a mass increase 
rate of 

— = 4.4- 11 X 10^12 yr-^ (3) 

In fact, in Ref. [32] it is argued that the DM content of the Sun cannot be 
due to a continuous accretion throughout its lifetime; however, the value of 
the amount of Galactic DM encountered by the Sun in the latest 4.5 Gyr in 
Ref. [12] radically differs from that in Ref. [43] . 

Accretion of DM, in the form of Mirror Matter [36], by the Sun during its 
lifetime is envisaged in Ref. [39j. Indeed, according to eq. (7) of Ref. ^J, re- 
produced below, a gravitating body of mass M traveling at speed v through 
a volumetric distribution of DM with density g will accrete an amount AM 
over a time span t 



M \MqJ VlO km s^V ylO"^^ g J \loio yj.^ 

Thus, if we assume v = 30 km s~^ and the upper limit of the local DM 
density in the solar system given by eq. the Sun will increase its mass 
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at a rate 

— ~ 7 X 10-12 yr-\ (5) 

close to eq. ([3]) from Ref. |42j. On the other hand, the mean Galactic value 
of eq. ([2]) applied to eq. (jH yields a mass accretion rate of 

— ~ 1.5 X IQ-^^ yT-\ 6) 

For several possible observable manifestations of Mirror Matter at various 
astronomical scales, see Ref. |47] . 

Concerning the planets, according to Table II and Table III of Ref. |43] . 
Jupiter would have captured the largest amount of DM during the latest 4.5 
Gyr followed by Saturn, Uranus and Neptune. See Table [1] for the related 
relative mass variation rates M/M of the eight planets. 



Table 1: DM accretion rates M/M on major bodies of the solar system 
according to Table II and Table III of Ref. ^43j . 

Planet M/M (yr^^) 



Mercury 


2.7 X 10- 


-17 


Venus 


1.6 X 10" 


-17 


Earth 


1.4 X 10" 


-17 


Mars 


4.1 X 10" 


-17 


Jupiter 


0.6 X 10" 


-17 


Saturn 


1.1 X 10" 


-17 


Uranus 


2.9 X 10" 


-17 


Neptune 


3.5 X 10" 


-17 



The estimates in Ref. [21] are different but, as we will show, such dis- 
crepancies will not substantially affect the main result of our analysis. 

We will investigate if the phenomenon of DM capture and accretion by 
the Sun and the planets has somewhat influenced the solar system dynamics 
throughout its 4.5 Gyr lifetime. We will address this problem in a specific 
way, i.e. we will look at the possible modifications of the orbital trajectories 
of the satellites of the major bodies which experienced the DM accretion, 
assumed not self-annihilating [45j. 

The paper is organized as follows. In Section [2] we will briefly review 
some basics of orbital motion around a mass- varying body. In Section [3] we 
will apply our results to concrete astronomical scenarios in the solar system. 
Section [His devoted to the Conclusions. 



4 



2 Orbital dynamics around a mass-varying body 



The acceleration experienced by a test particle orbiting a mass- varying cen- 
tral body with //(t) = GM(t) can be approximated with 



l+(^)(t-to) 



(7) 



where to is a given epoch and /i = /ul^^^g , A = Alt=to • GUI' case 
coincides with the present time, and t < Iq; thus, t — Iq < 0. Since in our 
case fi is due to DM accretion, we will assume that fi =cost. throughout 
the solar system's lifetime, i.e. we will not deal with /i; the validity of such 
an approximation will be justified a-posteriori in Section [3l Moreover, the 
condition of validity of eq. ([7|), i.e. — tg) !> is satisfied for all 

the major bodies of the solar system over At = —4.5 Gyr, given the DM 
accretion rates quoted in Section [TJ 

In Ref. [38] the secular, i.e. averaged over one orbital revolution, orbital 
effects induced by eq. ([7]) on the usual Keplerian orbital element^ of a test 
particle have been analytically worked out by means of the Gauss equations 
for the variations of elements. They are 



(1-e) V/^ 

(e) = {l + e)(^] (9) 



/) = (10) 



nj = (11) 

(dj) = (12) 
m) = n + 27r(^^y (13) 

where a,e,I,Q,u;,A4 are the values at the epoch of the semimajor axis, 
the eccentricity, the inclination, the longitude of the ascending node, the 
argument of pericenter and the mean anomaly, respectively; n = \J n/a? 
is the Keplerian mean motion. Such results are not in contrast with the 
orbital evolution effectively followed by the test particle, as shown in detail in 
Ref. [H]. Indeed, for, e.g., /i/// < 0, i.e. for a mass decrease, the trajectory, 
which is not an ellipse, would expand in such a way that the osculating 



^The type I, according to Ref. 
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semimajor axis and eccentricity of the Keplerian ellipses approximating it at 
the pericentre, which remains fixed, get reduced revolution after revolution, 
in accordance to eq. ([8]) and eq. ([9]) ; see Figure 1 of Ref . [38] and the related 
discussion. Moreover, also the osculating Keplerian period Pb = 27r/n, 
i.e. the time required to fully describe an osculating Keplerian ellipse, gets 
shorter, according to 

(^>^|-^^> 

while the time of two consecutive pericentre crossings along the true path 
grows. 

The variations of the radial p and transverse r components of the radius 
vector, taken at two consecutive pericentre passages, are [48] 

Ap = _^a(l-e)f^y (15) 



A. = '-^a(t)J\±I, (16) 

so that the overall pericenter-to-pericenter orbit variation is 



Ar = Ap,ll + A.^^±±. (17) 



3 Consequences of DM accretion on the orbits of 
the sateUites of the major planets 

Here we will fruitfully apply the results of Section [2] to the problem of de- 
termining the impact of DM accretion on the Sun and the major planets on 
the orbital evolution of their planetary and satellite systems throughout the 
past 4.5 Gyr. 



3.1 The Sun and the planets 

Concerning the Sun, by assuming for it M/M = 4.4 — 11 x 10^^^ yr~^ [l2], it 
turns out that 4.5 Gyr ago the orbits of the planets were larger than now by 
the amounts listed in Table [2j It can be noted that, while the distances of 
the inner planets at their formation are not substantially altered, the giant 
gaseous planets are shifted by about 1 — 10 au with respect to their present- 
day locations. If the rate of DM accretion on the Sun will be the same for the 
remaining lifetime of the solar system, the consequences for its planets would 
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Table 2: Increment Ar of the planetary orbits for M /M = 4.4 — 11 x 10 
yr~^ j42j and At = —4.5 Gyr according to eq. ()17p . 

Planet Ar (au) 



Mercury 


0.06 


- 0.15 


Venus 


0.09 


- 0.23 


Earth 


0.13 


-0.32 


Mars 


0.21 


-0.52 


Jupiter 


0.68 


- 1.71 


Saturn 


1.26 


- 3.17 


Uranus 


2.53 


- 6.31 


Neptune 


3.82 


- 9.55 



be relevant because of the continuous shrinking of their orbits. Indeed, when 
the Sun will have reached the tip of the Red Giant Branch (RGB) in the 
Hertzsprung-Russel (HR) diagram in the next 7.58 Gyr [50j, the orbits of the 
planets will be reduced by the amounts listed in Table El In particular, the 

Table 3: Decrement Ar of the planetary orbits for M/M = 4.4 — 11 x 10^^^ 
yr~^ [42] and At = +7.58 Gyr according to eq. (fTT]) . 



Planet 


Ar (au) 




Mercury 


-(0.10 - 


0.25) 


Venus 


-(0.15 - 


0.38) 


Earth 


-(0.21 - 


0.54) 


Mars 


-(0.35 - 


0.88) 


Jupiter 


-(1.15 - 


2.89) 


Saturn 


-(2.13 - 


5.33) 


Uranus 


-(4.25- 


10.65) 


Neptune 


-(6.43- 


16.08) 



heliocentric distance of the Earth may be reduced by ~ 0.2 — 0.5 au. Figure 
2 of Ref. [50] shows that at the beginning of RGB the solar photosphere will 
reach about 0.5 — 0.6 au; after entering the RG phase things will dramatically 
change because in only ~ 1 Myr the Sun will reach the tip of the RGB phase 
expanding up to 1.20 au, while the Earth's distance will be just 0.8 — 0.5 
au or less. However, it must be pointed out that the results of Ref. ^0] do 
not take into account possible DM accretion which, in principle, may have 
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consequences on the lifetime cycle of main sequence stars depending on the 
type of DM involved [51]. Another interesting consequence of the Sun's 
DM accretion is that the semimajor axis of the Earth's orbit, which 
can be approximately taken equal to the Astronomical Unit amounting tcH 
00 = 1.00000018 au, undergoes a secular increment of just 0.02 — 0.05 m 
yr^^ according to eq. dS]) and M /M = AA — 11 x 10^^^ yr"^- An anomalous 
increase of the Astronomical Unit of 0.15 ±0.04 m yr~^ was reported for the 
first time in Ref. ^52j. Later results are 0.07 ± 0.02 m yr~^ [53j and 0.05 m 
yr~^ [54j. Anyway, we point out that the figures presented here would be 
radically smaller if eq. ([2]) for the Galactic mean density was used in eq. (jH) 
of Ref. |39j . An independent check of the assumption on which the previously 
results are based, i.e. the Sun's mass accretion rate of eq. ([3]), comes from 
the upper bounds on the secular variation of the Newtonian constant of 
gravitation G which have been recently obtained by processing long data 
sets of planetary observations with the latest ephemerides. Actually, they 
have been derived from the variation of the gravitational parameter /i = GM 
of the Sun, but 

ii G M 

SO that they are useful for our purposes as well. W.M. Folkner, by analyzing 
the Mars spacecraft data set with the JPL DE421 ephemerides recently 
obtained [55] 

^ < 2 X 10"^^ yr"^ (19) 
Lunar Laser Ranging yields [56] 

= (4 ± 9) X 10"^^ yr"\ (20) 

with an upper bound of 1.3 x 10^^^ yr^^, and [S7] 

^ = (2±7) X 10-13 yr-\ (21) 

yielding an upper bound of 9 x 10^^'^ yr^^- Let us mention that E.V. Pitjeva, 
with the EPM ephemerides, recently obtained a statistically significative 
non-zero result [SH] 

= (-5.9 ± 4.4) X 10~i^ yr-^ (22) 



See http://ssd.jpl.nasa.gov/txt/p_elem_t2.txtlon the WEB. 
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the upper bound is —1.5 x 10^ yr • Such results contradict the possibihty 
that a steady mass accretion as large as that by eq. Q-eq. ([5]) may occur 
for the Sun; on the contrary, the smaller value by eq. ([6]) due to the mean 
Galactic DM density would be allowed. 

3.2 The planets and their satelhtes 

Let us, now, move to the planetary systems of satellites. In the case of the 
Earth and Moon, eq. ([IT]) and M/M = 1.4 x 10"^'^ yr^^ from Ref. [13] yields 
an expansion of the lunar orbit of just 0.16 km. The recently observed secular 
increase of the eccentricity of the geocentric Moon's orbit amounts to (e) = 
(0.9 zb 0.3) X 10~^^ yr^^ [59j. It is several orders of magnitude larger that 
the secular increase predicted by eq. ([9]). Concerning the observationally 
determined secular decrement a = —0.4 m yr~^ of the semimajor axis of 
the orbit of the geodetic satellite LAGEOS [6D], DM accretion is unable to 
explain it because eq. ([8]) yields an increase of the order of lO"^'^ m yr~^. 

In Table|3]-Table[7]we quote the expansion of the orbits of some satellites 
of Jupiter, Saturn, Uranus and Neptune 4.5 Gyr ago induced by the DM 
accretion according to Ref. |43J. It can be shown that the effect of such 

Table 4: Jupiter: increment Ar of the Galilean satellites' orbits and of 
S/2003J2 for M/M = 0.6 x lO"^'^ yr"^ [43] and At = -4.5 Gyr according 
toeq. (fTTl). 



Satellite Ar (km) 



lo 


0.07 


Europa 


0.11 


Ganymede 


0.17 


Callisto 


0.31 


S/2003J2 


7.21 



form of M/M is very tiny, amounting approximately to 10~^ — 10^ km. 
This shows that the discrepancies between the evaluations of Ref. |43j and 
Ref. |21] are negligible in our case. Moreover, a-posteriori we can well justify 
the assumption of neglecting /i. Let us note that DM accretion on Jupiter 
cannot be the cause of the observationally determined acceleration of the 
orbits of lo, Europa and Ganymede [61j, because the orbital period of lo is 
shortening, while those of Europa and Ganymede are lengthening; eq. (|14p. 
instead, predicts an uniform increase for Pb- We can conclude that DM 



9 



Table 5: Saturn: increment Ar of some satellites'orbits for M /M = 1.1 x 
10~^^ yr~i [43j and At = -4.5 Gyr according to eq. (fTTl) . 



Satellite 


Ar (km) 


Mimas 


0.06 


Titan 


0.39 


lapetus 


1.15 


Phoebe 


4.78 


Fornjot 


9.73 



Table 6: Uranus: increment Ar of some satellites'orbits for M /M = 2.9 x 
10~^^ yr~^ [33j and At = —4.5 Gyr according to eq. (fT7|) . 

Satellite Ar (km) 

Ariel 0.15 

Oberon 0.48 

Margaret 26.4 

Setebos 28.18 

Ferdinand 25.62 



Table 7: Neptune: increment Ar of some satellites'orbits for M /M = 3.5 x 

10-1^ yr-i [33] and At = -4.5 Gyr according to eq. (fTTl) . 

Satellite Ar (km) 

Naiad 0.05 

Galatea 0.06 

Triton 0.35 

Nereid 14.46 



accretion on the major planets of the solar system had negligible impact in 
sculpting the orbital patterns of their satellite systems. 

4 Conclusions 

In this paper we have investigated the impact that the phenomenon of DM 
capture by the Sun and the major bodies of the solar system may have on the 
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orbital dynamics of the planets and their satellites on timescales comparable 
with the age of the solar system. 

Concerning the planets, a solar mass increase by DM-assumcd not-sclf- 
annihilating-as large as « 10~^^ yr"^ would displace their orbits outward 
by about 10~^ — 10^ au at the epoch of the solar system's formation 4.5 
Gyr ago; for the Earth the increase of its heliocentric distance would be 
in the range 0.1 — 0.3 au, while the largest shifts are for the giant gaseous 
planets. The solar DM accretion rate is assumed by postulating that the 
maximum amount of DM in the Sun, evaluated as 2 — 5% of the total solar 
mass, comes from a continuous accretion process. Anyway, we note that 
the real occurrence of such a process is controversial; moreover, also in this 
case, it would be orders of magnitude smaller if the local DM density was 
as large as the mean Galactic one. On one hand, one of the consequences of 
the assumed DM-induccd solar accretion is that the osculating semimajor 
axis of the Earth's shrinking orbit would experience a secular increase as 
large as the latest observational determinations of the Astronomical Unit 
rate of 0.07 — 0.05 m yr""^. In fact, no real contradiction exists because it 
has been shown that the osculating semimajor axis of a gradually shrinking 
trajectory around a mass-increasing central body gets, in fact, larger. By 
assuming that the solar DM accretion will take place in the next billions 
year at the same rate as in the past, the orbits of the planets will shrink by 
10~^ — 10^ au (0.2 — 0.5 au for the Earth), with consequences for their fate, 
especially for the inner planets. On the other hand, long observational data 
records processed with the latest lunar and planetary ephemerides set upper 
bounds on possible secular variations of the Sun's gravitational parameter 
GM which arc more than one order of magnitude smaller than 10~^^ yr"^- 

The DM accretion on the planets has, instead, no effects on the dy- 
namics of their satellites whose planetocentric distances are altered by just 
10~^ — lO"*^ km over 4.5 Gyr. Concerning the Earth, both the observationally 
determined secular decrease of the semimajor axis of the LAGEOS satellite 
of 0.4 m yr~^ and the secular increase of the eccentricity of (0.9±0.3) x 10'^^ 
yr~^ of the lunar orbit cannot be explained by the DM accretion on the 
Earth. The same occurs for the observed accelerations of the orbits of lo, 
Europa and Ganymede. 
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